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Abstract 

On the basis of the gauge principle of field theory, a new variational formulation 
is presented for flows of an ideal fluid. The fluid is defined thermodynamically by 
mass density and entropy density, and its flow fields are characterized by symmetries 
of translation and rotation. The rotational transformations are regarded as gauge 
transformations as well as the translational ones. In addition to the Lagrangians rep- 
resenting the translation symmetry, a structure of rotation symmetry is equipped 
with a Lagrangian Ka including the vorticity and a vector potential bilinearly. Eu- 
ler's equation of motion is derived from variations according to the action principle. 
In addition, the equations of continuity and entropy are derived from the varia- 
tions. Equations of conserved currents are deduced as the Noether theorem in the 
space of Lagrangian coordinate a. Without A^, the action principle results in the 
Clebsch solution with vanishing helicity. The Lagrangian yields non-vanishing 
vorticity and provides a source term of non- vanishing helicity. The vorticity equa- 
tion is derived as an equation of the gauge field, and the A.a characterizes topology 
of the field. The present formulation is comprehensive and provides a consistent 
basis for a unique transformation between the Lagrangian a space and the Eulerian 
X space. In contrast, with translation symmetry alone, there is an arbitrariness in 
the transformation between these spaces. 

Key words: Gauge priciple. Variational formulation. Ideal fluid, Rotation 
symmetry, Chern-Simons term 
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1 Introduction 



A guiding principle in physics is that physical laws should be expressed in a 
form that is independent of any particular coordinate system. Fluid mechanics 
is a field theory of Newtonian mechanics that has Galilean symmetry: that is to 
say, it is covariant under transformations of the Galilei group. Two symmetries 
{i.e. transformation invariances) are known as subgroups of the Galilei group: 
translation (space and time) and space-rotation. In the present formulation of 
flows of an ideal fluid, we seek a scenario which has a formal equivalence with 
the gauge theory in physics. The gauge theory provides a basis for reflection 
on the similarity between fluid mechanics and other physical fields. 

The gauge principle ^ requires a physical system under investigation to have a 
symmetry, i.e. a gauge invariance (an invariance with respect to a certain group 
of transformations). Following this principle, an attempt has been made in 
Kambe (2007) to study a gauge symmetry of flow fields with respect to trans- 
lational transformations. The formulation started from a Galilei-invariant La- 
grangian of a system of point masses which is known to have global gauge sym- 
metries with respect to both translation and rotation. ^ It was then extended 
to flows of a fluid, a continuous material characterized with mass density and 
entropy density. In addition to the global symmetry, local gauge invariance 
of a Lagrangian is required for such a continuous field. Thus, the convective 
derivative of fluid mechanics {i.e. the Lagrange derivative) is identified as a 
covariant derivative, which is a building block in the framework of gauge the- 
ory. Based on this, appropriate Lagrangian functionals are defined for motion 
of an ideal fluid. Euler's equation of motion is derived from the action princi- 
ple. It is found that a general solution of this formulation is equivalent to the 
classical Clebsch solution (Clebsch 1859; Lamb 1932), in which the vorticity 
has a particular form and the helicity vanishes. 

In most traditional formulations, the continuity equation and isentropic con- 
dition are taken into account as constraints for the variations of the action 
principle by using Lagrange multipliers, where the isentropy means that each 
fluid particle of an ideal fluid keeps its entropy value along its trajectory (but 
the fluid is not necessarily homentropic^ .) In the present formulation however, 
those equations are derived from the action principle. The total Lagrangian 
density consists of kinetic energy and internal energy (with negative sign), sup- 
plemented with additional two terms associated with translation invariances 

^ See Appendix, or standard texts such as Weinberg (1995), Frankel (1997) or 
Aitchison & Hey (1982). 

^ These symmetries are interpreted as the homegeneity and isotropy of space in 
§7 and 9 of Landau & Lifshitz (1976). See §2.1 and Appendix for the definition of 
global gauge symmetry {i.e. global gauge invariance) and local gauge symmetry. 
^ A fluid is homentropic if its entropy per unit mass is uniform throughout space. 



2 



of mass and entropy (see (27) of §3.1). As is interpreted in detail in §4.1, it 
is remarkable that the Lagrangians associated with the latter two terms are 
given such forms as can be eliminated from the action represented with the 
particle coordinate a {i.e. Lagrangian coordinate), so that the Euler-Lagrange 
equations represented by the independent varibales a = {a^,a'^,a^) are unin- 
fluenced by the additional terms. Such a form of Lagrangian has close analogy 
to the Chern-Simons term in the gauge theory. 

Earlier papers (Kambe, 2003a, b) investigated the rotation symmetry of the 
velocity field v{x) in a preliminary way, and found that the vorticity u: = Vxv 
is in fact a gauge field associated with the rotation symmetry. However in the 
previous studies including Kambe (2007), the vorticity and helicity were not 
taken into account in satisfactory manner. Present formulation tries to account 
for the vorticity and helicity more properly and comprehensively. Thus, the 
main theme of the present paper is the rotation symmetry. 

A new structure of the rotation symmetry is given in §6 (after a review of the 
rotational transformations in §5) by the following Lagrangian A^: 



where A is a vector potential and = dtOJ + {v-'V)uj — {uJ-'V)v + (V-v)(jj. 

This is derived from a Lagrangian form characteristic of a topological term 
known in the gauge theory. This new term yields non-vanishing rotational 
component of the velocity field, and provides a source term of non-vanishing 
helicity. This is also the Chern-Simons term, describing non-trivial topology 
of vorticity field, i.e. mutual linking of vorticity lines. The vorticity equa- 
tion is derived as an equation for the gauge field. These are described in §6 
(Lagrangian associated with rotation symmetry), §7 (Corrections to the 
Variation), and §8 (Significance of A^). 

With regard to the variational formulation of fluid flows, the works of Eckart 
(1938, 1960) and Herivel (1955) are among the earliest to have influenced 
current formulations. Their variations are carried out in two ways: i.e. a La- 
grangian approach and an Eulerian approach. In both approaches, the equa- 
tion of continuity and the condition of is entropy are added as constraint con- 
ditions on the variations by means of Lagrange multipliers. The Lagrangian 
approach is also taken by Soper (1976). In this relativistic formulation those 
equation are derived from the equations of current conservation. Several action 
principles to describe relativistic fluid dynamics have appeared in the past (see 
Soper (1976, §4.2) for some list of them), not only that of Taub (1954). 

Firstly, in the Lagrangian approach, the Euler-Lagrange equation results in 
an equation equivalent to Euler's equation of motion in which the acceleration 
term is represented as the second time derivative of position coordinates of 
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the Lagrangian representation. In this formulation, however, there is a certain 
degree of freedom in the relation between the Lagrangian particle coordinates 
(a^,a^,a^) and Eulerian space coordinates {x^,x^,x^). Namely, the relation 
between them is determined only up to an unknown rotation. This point has 
not been considered seriously in the past, although it is fundamental for the 
formualtion of flows. This will be discussed more precisely in the sections 
4.2 and 8.2. A common feature of Herivel (1955) and Eckart (1960) is that 
both arrive at the Clebsch solution in the Eulerian representaion (see Sec. A 
and 3.2). However, it was remarked by Bretherton (1970) that the Clebsch 
representation has only local validity in the neighbourhood of a chosen point 
if vortex lines are knotted or linked. 

In the second approach, referred to as the Eulerian description, the action 
principle of an ideal fluid results only in potential flows if the fluid is homen- 
tropic, as shown in Kambe (2003a, 2007). However it is well-known that, even 
in such a homentropic fluid, it should be possible to have rotational flows. ^ 
This is a long-standing problem (Serrin 1959; Lin 1963; Seliger & Whitham 
1968; Bretherton 1970; Salmon 1988). Lin (1963) tried to resolve this dif- 
ficulty by introducing a constraint as a side condition, imposing invariance 
of Lagrangian particle labels (a^,a^,a^) along particle trajectories. The con- 
straints for the variation are defined by using Lagrange multipliers (functions 
of positions), which are called 'potentials'. In addition, the continuity equa- 
tion and isentropic condition are also taken into account by using Lagrange 
multipliers. However, the physical significance of those potentials introduced 
as the Lagrange multipliers is not clear. Mysteriously, the Lagrange multiplier 
for the continuity equation becomes the velocity potential for the irrotational 
part of the velocity. Gauge theory for fluid flows provides a crucial key to 
resolve these issues, which is the main target of the present investigation. 

Symmetries imply conservation laws. According to the relativistic formulation 
{e.g. Soper 1976; Landau & Lifshitz 1987 (Ch.XV)), as long as the flow fields 
obey the Euler-Lagrange equation, an energy-momentum tensor T'^" must 
obey the conservation equation of the form d^, T^^ = (considered in Sec. 4.2 
below). The Energy- momentum tensor T'^" in the present case is given in 
Sec. 4. A similar field-theoretic approach is taken in Jackiw (2002) by applying 
the ideas of particle physics to fluid mechanics in terms of Hamiltonian and 
Poisson brackets, both relativistically and nonrelativistically, and extension 
to supersymmetry is also considered. In Soper (1976) and Jackiw (2002), the 
nonrelativistic part follows the traditional approach and no explicit gauge- 
theoretic consideration is given to fluid mechanics in the sense of the present 



^ In fact, Euler (1755) showed that his equation of motion can drescribe rotational 
flows, in the historical paper: L. Euler, Principes generaux du mouvement des fluides, 
MASB, V. 11 (1757), 274-315. (MASB: Memoires de I'Academie Ryale des Sciences, 
annee 1755. Berlin 1757.) 
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formulation. 



The appendix describes the background of the present theory: (A) Known 
gauge invariances in two systems of electromagnetism and Clebsch form of 
fluid flows, (B) Related aspects in quantum mechanics and relativistic theory, 
and (C) Brief scenario of gauge principle. 



2 Translation symmetry of flow fields 

The formulation of Kambe (2007) started from a Galilei-invariant Lagrangian 
of a system of point masses, which is known to have global gauge symmetries 
with respect to both translation and rotation. Concept of local transformation 
is an extention of the global transformation to describe a continuous field. In 
order to represent a continuous distribution of mass, the discrete positions 
of point masses are replaced by continuous parameters a = (a^,a^,a^) of 
particles in a bounded sub-space M (of consideration) in the three-dimensional 
Euclidean space E^. The spatial position of each mass particle labelled by a 
(Lagrange parameter) is denoted by a; = Xa(t) = X{a,t), a function of a as 
well as the time t. An Eulerian space point is denoted by 3). The 

fluid particle occupying the point x at time t is denoted by a{x, t). ^ 



2.1 Global invariance of a fluid 

For a continuous distribution of mass {i.e. for a fluid), Lagrangian functional 
is defined as 



where x and v denote Xj^it) and v'^ = dtX^(t) {k = 1,2,3) respectively, and 
L{x,v) is a Lagrangian density. Suppose that an infinitesimal transformation 
is expressed by x ^ x' = x + 6x and v ^ v' = v + 6v, where 6x = ^{x, t) and 
Sv = dt{Sx), a component of ^ is denoted by ^*^(x, t), an arbitrary differentiable 
variation field. The resulting variation of the Lagrangian density L[x, v) is 




^ Here the Lagrangian coordinates a = (a^,a^,a^) are defined so as to represent 
the mass coordinate such that d^a = pd^x, where d^a = da^da^da'^ and d^a; = 
dx^dx^dx^ . A variable such as Xa denotes the one with respect to the particle a. 
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This does not vanish in general. In fact, assuming the Euler-Lagrange equation, 
dLjdx^ — dt{dL/ dv^) = 0, we obtain 




Global invariance, i.e. 6A = J 6L(x,v) d^x = for arbitrary constants ^'^ 
{k = 1,2,3), requires that the total momentum defined by j{dL/dv^)d?x 
must be invariant. In the local transformation, however, the variation field 
depends on the time t and space point x, and the variation 5A = / 5Ld?x 
does not vanish in general. 

2.2 Covariant derivative 

According to the gauge principle of Appendix C, non-vanishing of 6A is un- 
derstood as meaning that a new field G must be taken into account in order 
to achieve local gauge invariance of the Lagrangian A. To that end, we try to 
replace the partial time derivative dt by a covariant derivative D^, where the 
derivative Dj is defined by := dt + G. In dynamical systems like the present 
one, the time derivative is the primary object to be considered in the analysis 
of local gauge transformation, as implied by the expression (2). Thus, the time 
derivatives dt^ and dtq are replaced by Dj^ = dti + G ^ and 'Dtl = dtq + Gq, 
where q is understood to denote x. A free-field Lagrangian Af is defined by 

Af = 1 Li{v) d'x := \ J {v, v) pd?x = \j (D^g, D^g) d^a, (3) 

where pd?x = d^a denotes the mass element in a volume element d^x of the 
x-space with p the mass-density, and v = Dtq. 

2.3 Galilean transformation 
2.3.1 Global transformation 

Symmetries to be investigated are the translation symmetry and rotation sym- 
metry, which exist in the discrete system. The Lagrangian Af defined by (3) 
for a continuous field has such global symmetries, inherited from the discrete 
system. It is a primary concern here to investigate whether the system of fluid 
flows satisfies local invariance. We consider the translation symmetry first, and 
in the second half of the present analysis we investigate the rotation symmetry. 

A translational transformation from one frame F to another F' moving with 
a relative velocity L/ is a Galilean transformation in Newtonian mechanics. 
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The transformation is defined by 



X 



(t, x) =^ x' = {t\ x') = (t, x-Ut). (4) 



Tliis is a sequence of global translations with a parameter t. Corresponding 
transformation of velocity is as follows: v = (1, v) =^ v' = (1, v') = 
(1, V — U). From (4), differential operators dt = d/dt and dk = d/dxk {d'^ = 
d/dx'i,) are transformed according to dt = dt' — U ■ V' and V = V', where 
V = {di, 82, ds) and V' = {d[, 82, d'^). There is a transformation invariance of 
the convective derivative Dc = dt + {v- V): namely dt + {v- V) = df + {v' ■ V'). 



2.3.2 Local Galilean transformation 

The velocity field v[x,t) is defined hj the particle velocity, i.e. 

v{Xa,t) = ^^Xa{t). (5) 

Consider the following infinitesimal transformation: 

x'{x,t) = x + ^{x,t), t' = t. (6) 

This may be regarded as a local gauge transformation between non-inertial 
frames. In fact, the transformations (6) is understood to mean that the co- 
ordinate a; of a fluid particle at a; = Xa{t) in the frame F is transformed 
to the new coordinate x' of F', given by x' = X'^{Xa,t) = Xa{t) + ^{Xa,t). 
Therefore, its velocity v = {d/dt)Xa{t) is transformed by 



v'{x')^^X', = + = viX,) + (d/dt)^„ (7) 



dt " dt 
{d/dt)^, = dt^ + {vV)^ 



X—Xa 



This is interpreted as follows: the origin of a local frame is displaced by — ^ 
where its coordinate axes are moving (without rotation) with the velocity 
— (d/dt)^ in accelerating motion (non-inertial frame). This implies that the 
velocity v{x) is transformed according to (7), since the local frame is moving 
with the velocity — (d/dt)^. Note that, in this transformation, the points x 
and x' are the same with respect to the space F. In view of the transformation 
(6), the time derivative and spatial derivatives are transformed as 



dt = d, + {dt^)-V\ V' = (9;), (9) 
dk = d',+ d,^,d\, d', = d/dxi. (10) 
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2-4 Invariance of Dt and definition of velocity 

In this local transformation too, there is an invariance of the convective deriva- 
tive Dc = dt + {v ■ V) = df + (v' ■ V'). In fact, from (10) and (7), we obtain 

. V = ■ V + ■ VO ■ V = v'{x') ■ V + ( - (d^/dt) + v-V^)- V, 

where v = v' — d^/dt. The last term is —dt$, ■ V by (8). Thus, we find the 
invariance of by (9). This invariance {i.e. the transformation symmetry) 
implies that the covariant derivative is in fact given by the convective 
derivative D^. 

Bt = dt + iv-V)= dt, + {v' ■ V) (11) 
(see Kambe (2007) for its direct derivation). 

For a scalar function f{x, t), we have Df/ = dtf + v ■ V/ (Lagrange derivative 
oi f{x, t)). It is important to recognize that this derivative Dtf has an intrinsic 
meaning. That is, we have the invariance: Dtf = D[,f', with respect to the 
local Galilean transformation (6), since a scalar means /' = f'{x', t') = f{x,t). 
In general, vectors or tensors need not have this property. 

By using the covariant derivative Dt, we can define the velocity v by DtX. 
In fact, each component a* of the particle label a is a scalar, and satisfies 
Dja* = dta"^ + (v ■ V)a* = 0, since the particle with the label a* moves with 
the velocity v = dtX{a,t) by the definition (5). Setting x = X{a,t) for the 
particle position, we have 

DtX{a, t) = dtX{a, t) + Dta ■ V aX = dtX{a, t) = v, (12) 

where (VaX) = {dX'' /da''). In addition, applying Df to (6) and using D^, = 
Dt, we have v' = D[,x' = Dt{x + £) = v + Dt^ (consistent with (7)). Thus, 
the particle velocity v can be defined by v{x,t) = DtX. 



2.5 Transformation of acceleration 

The gauge transformation of Sec. 2. 3. 2 defines the transformations: x ^ x' = 
X + and v{x) — > v'{x') = v{x) + D^^. Namely, the variations of and 
are given by 

6x' = ^Tk x' = e, {Tk = d/dx^) (13) 
5v' = DtC = dtC + v''m\ (14) 
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where T^x* = Si- Next, let us consider acceleration of a fluid particle. Differ- 
entiating (7) with respect to t again, we obtain 

|X = ^^.W + |«(X.,*). (15) 

This describes a relation between accelerations in the two frames F and F'. Let 
us deflne the covariant derivative of the velocity v = (f *) by Dtv = dtV+v^Tk v. 
Variation of dtV is given by 

5{dtv) = dt{5v) + {5dt)v = dt{Y)ti) - {dteTk)v , 

where 5dt = dt' — dt = ~dt^^ from (9) (to the flrst order of ^). Then, using 
5v = DtC, and 6Tk = T^—Tk = —{TkS})Ti from (10), the transformation of D^f 
is given by 6{Dtv) = S{dtv) + {5v^)TkV + v^{5Tk)v + v^Tk{5v), namely. 

Thus, we obtain 

5(DitO = Di(DiO, 'Dt = dt + v^Tu. 

It is found that the covariant derivative DtV transforms just like the accelera- 
tion A = (d/dtyXait) according to (15). This implies that D^f represents the 
particle acceleration A{x,t): 

A = Btv = {dt + v''dk)v. (16) 



2.6 Galilean invariant Lagrangian 



According to the formulation of Kambe (2007) based on the gauge principle, 
the total Lagrangian is deflned by 

At = Af + Ae + + A^, (17) 
as far as the translation symmetry is concerned, where 



Ai = J l{v{x),v{x)) pd^x, K = - j^^e{p,s) pd?x, (18) 
A^ = - f il)t<f))pd^x, K = -f {Bti^)psd'x. (19) 

The variables v, p, e and s are the velocity vector, density, internal energy (per 
unit mass) and entropy (per unit mass) of the fluid, (j){x,t) and ip{x,t) are 
gauge potentials associated with the mass density p and entropy density ps 
respectively. The velocity is deflned by v{x, t) = Djcc, where Dt = dt + v ■ V . 
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From §2.4, it is almost obvious that the two Lagrangians and of (19) are 
invariant with respect to the local Galilean transformation {i.e. local gauge 
transformation) . 



3 Variational principle: Translation symmetry 

3. 1 Action principle 

The action / and the action principle 5/ = are defined by 

SI = 6 r f dtd^x Lt = 0, (21) 

Jto JM 

Lt: = Ip{v,v) - p e(p, s) - p Dt0 - ps Dtip . (22) 

where Lt is the Lagrangian density. There is a certain thermodynamic prop- 
erty that must be taken into consideration. That is the definition of an ideal 
fluid. Namely, there is no dissipation of kinetic energy into heat, i.e. there is 
no heat production within the ideal fluid. By thermodynamics, if there is no 
heat production, we have T6s = (T: temperature). Then, 




where p is the fluid pressure, h is the enthalpy defined by e + p/p, and ( ■ )s 
denotes s being fixed. However, the entropy s may not be uniform and depend 
on X through the initial condition. 

One of the aims of the present analysis is to show that the previous formu- 
lation (Kambe 2007) based on the symmetry of parallel translation alone is 
not sufficient to represent rotational flows with non-vanishing helicity. The 
rotaional symmetry will be considered from §5 below. 

3.2 Outcomes of variations 

Writing the Lagrangian density (22) as Lt[v, p, s, 0, ■i/'], we take variations of 
the variables v, p, s and potentials and ip, where variations are assumed 
independent. Substituting the varied variables v + 6v, p + Sp s + 6s, + 6(f) 
and ip + 6ip into Lt:[v, p, s, 0, ip] and writing its variation as 6Lt, we obtain 
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SLt = 5v ■ p{v — Vcj) — s Vip ) — Ss p Dtip 
+ 6p ilv^-h-Dt^-sDt^) 

+ 6^(dtp+V ■ (pv)) -dt{p6<P)-V ■{pv8<P) 

+ 5ij(dt{ps) + SJ -{psv)) -dt{ps6^)-V-{psv6i>), (24) 



where {d/dp){pe) = e + p {de/dp)s = e + p/p = h is used. The variation fields 
are assumed to vanish on the boundary surface enclosing the domain M, as 
well as at both ends to and ti of time integration of the action I (where M is 
chosen arbitrarily), as is usually the case. Thus, the action principle SI = 
for independent variations 6v, 6p and 6s results in 

Sv: V = V (/) + sViJ , (25) 
Sp: Iv"^ - h-Dt(l) - sBtiJ = 0, (26) 
6s: Dtip = dtip + v-\/ip = 0. (27) 

Using (25) and (27), we have 

Dt(j) = dt(f) + v ■ V0 = dt(f) + V ■ {v — sVip) = f ^ + dt(j) + s dti>. 

Under this relation and (27), the equation (26) can be rewritten as 

+ h + dt(t) + sdti} = Q. (28) 

This is regarded as an integral of motion. From the variations of 5(j) and 5ip, 

5<P: dtp + V ■{pv)=Q, (29) 
5^: dt{ps) + \7 ■{psv)=Q. (30) 

Using (29), the second equation reduces to the adiahatic equation: 

9js + V ■ Vs = Dts = . (31) 

Thus, we have obtained the continuity equation (29) and entropy equation 
(31) from the action principle. 

It is shown in Kambe (2007) that the present solution is equivalent to the 
classical Clebsch solution (Clebsch 1859). With the velocity (25), the vorticity 
uj is given by 

oj = V X V = Vs X S7ip. (32) 

This implies that the vorticity is connected with non-uniformity of entropy. 
Using (28), we find that the following Euler's equation of motion is satisfied: 

dtv + LJ X V = -W (^v^ + h). (33) 
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In this case, the hehcity vanishes (Bretherton, 1970). For a general velocity 
field, the helicity H defined by 

H[V] = j u-vd^x (34) 

is a measure of linkage and knottedness of vortex lines, and does not vanish 
in general (see Sec.8.1). 

3.3 Homentropic fluid 

For a homentropic fiuid in which the entropy s is a uniform constant Sq at all 
points, we have e = e(p), de = {j>/ p^) dp and dh = (l/p)dp from (23) since 
6s = 0. In addition, the motion is irrotational. In fact, from (25), we have 
V = V $ where $ = + sq "i/^, and a; = from (32). The integral (28) becomes 
+ h + dt^ = 0. The Euler's equation (33) reduces to 

dtv + V{^v^) = -V h, where Vh=^Vp. (35) 

Note that the left hand side is the convective time derivative for the irrota- 
tional velocity v'' = dk^: 

dtV + \/{^v^) = dtV+{v-\/)v = DtV. (36) 

Thus, as far as the action principle is concerned for a homentropic fiuid, Euler's 
equation of motion reduces to that for potential fiows of a perfect fiuid. 

Traditionally, this property is considered to be a defect of the formulation of 
the Eulerian variation described in the previous section, because the action 
principle should yield the equations for rotational fiows as well. In order to 
remove this (apparent) fiaw, Lin (1963) introduced the condition for the con- 
servation of the identity of particles denoted by a = (a^), which is represented 
by an additional subsidiary Lagrangian of the form / Ak-DtO''' d^x. This intro- 
duces three potentials Ak{x,t) as a set of Lagrange multipliers of conditional 
variation. Problem is that physical significance of Ak is not clear (Bretherton 
1970; Seliger & Whitham 1968; Salmon 1988). 

Let us recall that we have been considering the Lagrangian At satisfying the 
symmetry of parallel translation. However, the fiow field has another symmetry 
of rotational invariance (Kambe 2003a, b). The equation (31) implies that the 
entropy s plays the role to identify each fiuid particle. Furthermore, local 
rotation is captured by the expression (32). However, in a homentropic fiuid, 
there is no such machinery to identify each fiuid particle. Gauge invariance 
with respect to local rotation could be a candidate instead of s. 
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Insofar as the flow field is characterized by the translation symmetry alone, 
we have arrived at the above result, i.e. the flow field should be irrotational ii 
the fluid is homentropic. The fluid motion is driven by the velocity potential 
where $ = + Sq "0 with (p and ijj the gauge potentials. 



4 Equations in a-space 



4^.1 Lagrangian 



Let us consider another variational formulation with the Lagrangian repre- 
sented by the particle coordinates a = (a^,a^,a^) = (a,6, c). Independent 
variables are denoted with a'^ [fi = 0, 1, 2, 3) with being the time variable 
written as r {= t). The letter r is used instead of t in combination with the 
particle coordinates {k = 1,2,3). The physical-space position of a particle 
is expressed by X'^i^a^), or (X,Y,Z). Its velocity v'^ is given by drX^ = X^. 
The derivative dr = d/dr is equivalent to the covariant derivative : 

dr = dt + udn, + vdy + wdz = dt + v - V . (37) 

Mass density p is defined by the relation, pd^x = d^a, where d'^x denotes 
the volume element dX^dX^dX^. Using the Jacobian determinant J of the 
transformation from a-space to X-space, we obtain 

J' d{a^,a'^,a^) d{a,b,c) 

The total Lagrangian (17) with (18) and (19) is rewritten as 



At = y i X^ X^ d^a - J e(p, s) d^a - J d^cj) d^a- J s d^^ d^a. (39) 
The action / is defined by integration of At with respect to the time r: 




In this definition, the third integral /s = / dr / dr'pdC'a can be integrated with 
respect to r, and expressed as /[0]d'^a, where [0] = 01^-2 — 0|ri is the difference 
of at the ends T2 and ri and is independent of r G [Ti,r2]. Likewise, the 
fourth integral can be expressed as /4 = j[ip\sd?a., because the entropy s is 
independent of r according to (31). This means that the gauge potentials 
and -0 do not appear in the variations of the action I for r G [ti, r2]. In other 
words, the gauge potentials and tfj make their appearance only in the action 
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represented with the physical space coordinates x. The fact that the last two 
Lagrangians in (39) do not contribute to the variations (for r G [ti, r2]) implies 
that local Galilean symmetry is in fact an internal symmetry. 

Thus, by omitting the integrated terms, we have 



I = I L{X^) d'a, L=\Xl X', - 6(Xf , a% (40) 

where = dX^/da\ X^ = X^ = v'' {k, / = 1, 2, 3), and d^a = dr d^a. In 
the a-space, the point a is motion-less, while in the a;-space, there is a flow 
of particles driven by the potentials and iIj according to (25). 



4-2 Noether's conserved currents 



The Euler-Lagrange equation can be writtten as 

{d^ = d/da^). The energy-momentum tensor is defined by 

where k = 1,2, 3. As long as (41) is valid together with drL = 0, we have a 
conservation equation di^Tj^ = 0. In fact, we have 



dL 



-d„x'^-d,{x'^) 



dX^ 



dXt 



■Xt 



dL 
dX^ 



This is the Noether theorem (Noether (1918), Weinberg (1995)). 

By using the Lagrangian L of (40), we obtain dL/dX^ = X^ = v'^. The 
internal energy e depends on the density p which in turn depends on Xf = 
dX^ /dd by (38), where the entropy s depends on a' only. Hence, we have 

dL de dp p dp 



dxt 



dp dX^ 



p2 dXl 



k ' 



(42) 
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by using the thermodynamic relation given in (23). In view of (38), we obtain 



dL 



dL 



dXl ^ 9(a2,a3) ' dXf ^' d{a?,a^) ' dXf ^ d{a\a^j 
Thus, the energy-momentum tensor are found as follows (Eckart 1960): 



d{X\X^) dL 



p'-m^^ (43) 



To° = X^^ j -L = X^^X"; -L = \v^ + e = H (energy density), (44) 



'~^\dxi 



9(r, 6, c) 



,2 „5(x,r,z) 



To=p- 



d{a, r, c) 



^0 =P 



{de/dp = p/p^] 

d{X,Y,Z) 
d{a, b, t) 



Introducing the symbols (a, /5, 7) to denote (a, b, c) cyclically, another twelve 
components are 



X^" — Xq,Xt- + YaY-r + ZaZr, 



T" = X, 



9^ \ j^_p d{X,Y,Z) 
J a(a,/5,7) 



X'^ 



° V9X^ 



1 „,2 



(45) 
(46) 



" ^ dX^ 



0, 



T2 = 0. 



The Noether conservation law d^T^^ = reduces to the momentum equations 
(see below) for n = a: 



drVa + 9a -F = (for a = a,b,c), 



(47) 



where F = —^v'^ + h, and other two equations are obtained with cyclic per- 
mutation of (a, fe, c). Integrating this with respect to r between the limits 
and t, we find the Weber's transformation (Lamb 1932, Art. 15): 



K(r) = X,X. + + Z^Z, = \4(0) - d^x, 



(4J 



where x = /o -^dr = Jo{—^v^ + h)dT. For /_i = 0, the conservation law 
d^T^^ = describes the energy equation: 



drH + da 



d{X,Y,Z) 
d{T, b, c) 



P- 



d{X,Y,Z) 
d{a, r, c) 



P- 



d{X,Y,Z) 
d{a,b,T) 



P- 



0. (49) 



In the Eulerian description with the independent variables {t, x, y, z) instead 
of (r, a, 6, c), this reduces to 



iyt{\v'' + h) = -dtp. 



1 

p 
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This can be transformed to the following form of conservation of energy: 



0. 



(50) 



Using (45) and (46), the equation (47) reduces to the equation for the accel- 
eration AJt, a): 



p 



(51) 



which is known as the Lagrangian form of the equation of motion (Lamb 1932, 
Art. 13). This is obtained by noting that dah = {l/p)daP and 

In view of (16), the equation (51) is equivalent to Euler's equation of motion 
in the Eulerian description of the barotropic relation h{p) = dp' / p{p'): 



dtv + (v ■ V)v = --Vp= -Vh. 
P 



(52) 



Velocities V^(t, a) in the a-space may be determined by (48) with given initial 
conditions of Va{0,a) and h{0,a) at a = (a, 6, c) = x. However, in order 
to transform it to the aj-space, there is some freedom. The equation (48) 
is invariant by rotational transformations of a displacement vector AX = 
{AX, AY, AZ) since the left hand side is of the form of a scalar product with 
respect to AX, so that the vector AX is not uniquely determined (the density 
p and enthalpy h are not changed by this transformation since it is volume- 
preserving). The same freedom applies to the acceleration Aa{T,a) of (51) 
as well. Rotational transformations are considered in the next section. This 
formulation helps to eliminate the arbitrariness. 



5 Rotational transformations 



Next, in order to equip the action / with an additional mathematical structure 
of rotational symmetry (to be considered in §6), we summarize the rotational 
transformations here, and consider the rotational gauge symmetry. The related 
gauge group is the rotation group 5*0(3). An infinitesimal rotation is described 
by the Lie algebra so(3) of three dimensions. The basis vectors of so(3) are 
denoted by (61,62,63), which satisfy the commutation relations: 

[ej, 6fc] = Bjkiei , (53) 
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where ejki is the completely skew-symmetric third-order tensor, and 










1 




0-10 


ei = 


0-1 


, 62 = 





, eg = 


1 




1 




-10 








A rotation operator is defined by ^ = O'^e^ where 6'' [k = 1, 2, 3) are infinites- 
imal parameters: 

-9^ 9^ 
9^ -9^ . 
-9^ 9^ 

Then, an infinitesimal rotation of the displacement vector s = is 
expressed by 9 s. This is also written as ^ x s, where 9 = {9^, 9"^, 9^) is an 
infinitesimal angle vector. 

Preliminary studies of the rotational symmetry of fiuid fiows are given in 
Kambe (2003a, b). However, more details are elaborated in the following sub- 
sections 5.1 ~ 5.3, and the section 6 proposes a new Lagrangian La to account 
for the rotation symmetry. 

5.1 Gauge transformation (Rotation symmetry) 

We consider local rotation of a fiuid element about an arbitrary reference point 
xq within a fiuid. The non-abelian property of the rotational transformation 
requires infinitesimal analysis. At a neighboring point Xq + s for a small s, we 
consider its local transformation, s — s> s + 6s. 

Rotational gauge transformation of a point s in a frame F (an inertial frame, 
say) to s' in a non-inertial frame F' is expressed as 

s'{s,t) = s + 5s{s,t), t' = t, (55) 

instead of (6) for the translation. The variation 6s is defined by an infinitesimal 
rotation 9 = 9{t) (skew-symmetric) applied to the displacement vector s: 
6s = r]{s, t) = 9s = 9 X s. The corresponding transformation of velocity 
{v = Dfs) is obtained from (55) as follows: 

v'{s') = Dts' = Dt{s + 6s) =v{s) + DtT], (56) 
6v{s) = Dt'n = dMs,t) + {v-Vs)r}{s,t) = dt9s + 9v, (57) 



9 = {9,j) = 9''ek = 
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instead of (7) and (8) where = (d/ds^). As before (Sec. 2.3.2), in this 
transformation the points s and s' are the same point with respect to the 
frame F. This means that the coordinate frame F' is rotated by an angle —6 
with respect to F with a fixed origin Xq. At the origin a^o (s = 0), we have 
the following: 

Sv(0) = ev = e''ekV, (58) 
from (56) and (57), where v = v[xo). 

5.2 Acceleration 

Let us consider the time-derivative dtV in the Euler equation (35) for a home- 
ntropic fluid. It is proposed that this term be replaced with d^w as 

dtV — > dtV = dtv + ilv, f2 f = il'^Ck V = tl X V, (59) 

where Q is an operator defined by f2 = Q'^Ck for scalar fields fl'^, and Cl is 
an axial vector defined by {Q^, fi"^).^ The tensor Qij = {Vt^ek)ij is skew- 
symmetric {Vtij = —Qji). 

We consider variations at the origin (s = 0). Then, we have 6v = O^CkV from 
(58). Variation of Vt^ is defined by^ 

Sn'' = dt.e'' + EkimO^ ■ (60) 

Using 6dt = df — dt = —dtrj^Tk from (9), we obtain the variation of dtv: 

5{dtv) = dt{6v) + {6dt)v + 6n'' ekv + Q'^Ck Sv = dt{e^ek v) - dt{d sfnv 

+ (9t^' + SkimO' n^)ek V + QJ^Cm e'ei v, 
where Ekim can be eliminated by using the following equality obtained from 
(53): Eimkek = eiCyn - CmCi. Setting s = 0, we obtain 

5{dtv) = e^Ck [dtv + VTcrn v] + 2dte^ CkV = e X dtv + 2dte X v. (61) 

The second term of the last expression denotes the Coriolis term in the rotating 
system. Thus, it is found that the covariant derivative dtV behaves like the 
acceleration in a rotating frame, as far as Q transforms according to (60). 

In the previous analysis for the translation symmetry, it was shown that the 
left hand side of (35), dtV + V{^v'^), is in fact equal to the particle acceleration 
DtV of (36). In the present section taking into account the rotation symmetry, 

^ A simple font v is used here for brevity instead of bold- face three- vector v. 
^ The vector form is dCl = dtO + 6 x (l. 
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the partial derivative dt is replaced with the covariant derivative = dt + 
Thus, the particle acceleration Vj'y should be defined by 

WtV = dtV + VLv + grad(i;V2) =dtV + VLxv + grad(wV2) , (62) 

where VL is the gauge field with respect to the rotation symmetry. 

5. 3 Vorticity as a gauge field 

It can be verified that the gauge field VL in fact coincides with the vorticity 

V X v by the requirement of Galilean invariance of the covariant derivative 

of (62), as follows. Under a global Galilean transformation from one frame 
F : X = {t, x) to another : = (t*, a;*) which is moving with a uniform 
velocity U relative to F, the position vector and velocity vector in the 
frame are given by Ut and = V — U . Transformation laws of 

derivatives are dt = dt, — U ■ V* and V = V,,. Applying these and replacing 

V by V* + U, the covariant derivative of (62) is transformed to 



= dt,v^ + X + V, {vl/2) - [U ■ V^)v^ + (lxU + V^{v^-U) , 

since is a constant vector and V*(f/^) = 0. We require the covariance, i.e. 
Vjf = (Vfv)*. Namely, the right hand side should be equal to dt,v^ + x 
V* + W^{yl/2). Therefore, 



Q = {(l-VL^)xv^-{U ■ V*)v* + (lxU + V^(v^-U) 
= {(l-n^)xv^ + {h-V^xv^)xU, (63) 

on using the following vector identity: U x{y ^xv^) = — ([/■ V*)i'*+V*(L/"-v*), 
for the constant vector U. Equation (63) is satisfied identically, if 

n = Cl^, n = V XV = V^xv^ = Cl^. (64) 

The second relation holds by the Galilean transformation since V = V* and 
V = v^: + U . Thus, the Galilean invariance of VtV results in 

n = V xv = uj. (65) 

Namely, the gauge field Cl coincides with the vorticity uj. Consequently, the 
covariant derivative V^v (particle acceleration) is given by 

VtV = dtv + V{lv^) + UJ X V. (66) 
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Using the vector identity 



V X {V xv) = V(||vn - {v ■ V)v 



(67) 



the last expression (66) is transformed to the convective derivative of v {i.e. 
the particle acceleration): 



Thus, we have Vt = dt + {v ■ V) . So far, xq + s was a local position vector 
around a fixed point Xq. In general, using x in place of a^o + we arrive at 
the same definition: v{x,t) = VtX = (dt + (v ■ V))£C. 



6 Lagrangian La associated with rotation symmetry 

Associated with the rotation symmetry, we try to introduce an additional 
Lagrangian according to the gauge principle. The property of I3 and 1^ con- 
sidered below the equation (39) suggests how to find such a Lagrangian. It 
was seen there that the Lagrangians represented with the particle coordinates 
a = (a^,a^,a^) = (a, 6, c) can be integrated with respect to the time r and 
eliminated in the action to deduce the Euler-Lagrange equation, while they 
are non-trivial in the physical-space coordinates x = (x^,x^,x^) = {x,y,z) 
because of the non-trivial Jacobian matrix [dx^ /do!') of the transformation. 
As a result, the action principle yields the continuity equation and entropy 
equation. Thus the gauge potentials are regarded as mathematical agents with 
mechanical functions that drive fluid motions in the physical a;-space, whereas 
they disappear in motionless a-space. 

It is important to observe that, in the Lagrangian (39), the integrands of the 
last two integarals are of the form dr{-), because the entropy s and mass 
element pd^a; = d^a are independent of r. The action is defined by / = 
/ j[L' + dr{-)] drd'^a, where L' denotes the part of the first two terms of (39). 
This property is the simplest case of representation of topology in the gauge 
theory. In the context of rotational flows, it is known that the helicity (or 
Hopf invariant, Arnold & Khesin (1998)) defined by (34) describes non-trivial 
topology of vorticity field, i. e. mutual linking of vorticity lines. This is closely 
related with the Chern-Simons term in the gauge theory. The Chern-Simons 
term lives in one dimension lower than the original four-space-time {x^) of the 
action / because a topological term in the action is expressed in a form of total 
divergence {d^F^^) and characterizes topologically non-trivial structures of the 
gauge field (Chern 1979; Jackiw 1985; Desser et al. 1982). Here, we follow the 
formulation of Sec. 4.1 and look for a r-independent field directly. 



VfV = dfV + (v ■ V) V. 



(68) 
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6.1 Lagrangian in a-space 



The r-independent field can be found immediately from Eq. (47). Taking the 
curl of this equation with respect to the coordinates (a, b, c), we obtain 

V„ X d^Va = driVa X Va) = 0, (69) 

since x VaF = 0, where Va = (<9a, db, d^), and Va = (K, K) is defined 
by (45). It is found that Va x Va is independent of r. ® Hence, one may write as 
Va X Va = Oa{a), where the right hand side is a vector depending on a only. 
It is useful to formulate this property in the framework of exterior differential 
forms and Lie derivative. Vector Va is understood as a transformed form of 
the velocity v = (X^, Yr, Zr) = {u, v, w) into the a-space. This is seen on the 
basis of a 1-form V defined by 

V = Vada + Vbdb + dc (written as Va ■ da) (70) 
= udx + V dy + w dz. (71) 

where Va = uXa + vi/a + wza, Xa = OX/ da, u = Xt-, etc. and da = (da, db, dc). 
Hence Va is equivalent to (45). Its differential dV gives a two-form = dV: 

fl'^ = dV = Qa d6 A dc + Qb dc A da + da A db = fta ■ 

= dy A dz + Uy dz A dx + uj;^ dx A dy = u! ■ s^, (72) 

where Va x Va = {fla,^b,^c) = ^a, and V x t; = (ujx,ujy,ujz) = (j^J is the 
vorticity. Thus, it is seen that fia is the vorticity transformed to the a-space, 
and the 2- forms and S'^ are surface forms defined in the footnote ^ . The 
equation (69) is transformed to the r-derivative of the 2-form Q^, which is 
understood as the Lie derivative: Cg^Q"^ = 0. Its representation in the {x,t) 
space is given as follows. 

It is useful to define two tangent vectors W and V (tangent to the flow gener- 
ated hj v): 

W = dt + V, V = ud^ + vdy + wd, = v'^dk = vV, (73) 

where V is regarded as a tangent vector dual to the 1-form V of (71). The Lie 
derivative Cg^ can be defined by 

^dr = = ^dt+V = dt + Cy, (74) 

® This property can be related to the invariance of the Lagrangian with respect to 
the relabeling transformation a a', satisfying the mass invariance d'^a' = d^a. 
^ = (dx^Adx^, dx^Adx\ dx^Adx^), and 5^ = (da^Ada^, da^Ada^, da^Ada^). 
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according to the symbols of differential calculus and tangent vector, where 
V = v^dk is also a differential operator along a streamline generated by the 
flow {u,v,w). By applying the Lie derivative Cq^ on the 2-form Q"^ (Frankel 
1997, Ch.4.2), we obtain 

= Cd^n^ = dtn^ + Cv^^ = dtOJ ■ + [v, uj]-s^ + {y ■ v)u ■ 

= Cw^^ =Ew[i^] ■ s^ (75) 

where [v, oj] is the Lie bracket defined by (v ■ V)a; — (uj ■ V)v, and 

^,^[0;] = dfOJ + [v, uj] + (y ■ v)u: = dtuj + V x {u: x v), (76) 

where V ■ = is used. Thus, we obtain the vorticity equation from (75): 

dtu: + V X {u: X v) = 0. (77) 

Next, let us introduce a gauge vector-potential Aa = {Aa, Af,, Ac) in the a- 
space, and define its 1-form A^ by 

A^ = Aada + Ab db + Ac dc (written as Aa ■ da) (78) 
= Ax dx + Aydy + A^ dz. 

The exterior product of A^ and fi^ results in a volume form d^a = da A {db A 
dc) = db A {dc A da) = dc A {da A db) multiplied by a scalar {Aa, fla)'- 

A^An^ = {Aa ■ da) A {Qa ■ S^) = {Aa, Qa) d=^a , (79) 

where {Aa, fla) = Aa^la + AbQb + Ac^lc is a scalar product of Aa and Qa- 

In the X = {x, y, z) space, the same exterior product is given by 

A^ An'^ = {A, uj) d^x , d^x = dx A dy A dz. (80) 

It is obvious that the scalar product {A, uj ) is invariant under local rotational 
transformations in the x space. 

Thus, it is proposed that a possible type of Lagrangian is of the form, 

{Aa, na ) d^a] , 

where the time derivative dr was taken out of the inner-product since fla is 
independent of r. It follows that can be dropped in the action / expressed 
with the coordinate a except the integrated terms depending on the initial 
and final times of r-integration. However, this Lagrangian becomes non-trivial 
when transformed to the physical-space coordinate x because of the non-trivial 
Jacobian of transformation. It is shown later that Aa may be understood as 
a kind of vector potential associated with the gauge field fla- 



Aa = - {drAa, na)d'a=- Lq^ 
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6.2 Eulerian coordinates 



From the equality of (79) and (80), we have (A^, ria)d^a = {A, a;)d'^a;. Using 
the definition (74), 



dr 



w 



Cw 



where (75) was used. With the same calculus as in (80), the right hand side is 
given by [Cw A = {CwA, uj) dc'x, where the Lie derivative of a 1-form 
A^ is defined by 



CwA^ := {dtA, + v'^dkA, + AAv^) dx' = dx\ (81) 
:= {CwA), = dtA, + v^duA, + AAv'' (82) 

{e.g. Frankel (1997), §4. 2c) . Thus, the Lagrangian Ka in the x variables is 
Aa = -[ {CwA,u)d^x = -( a;)d^a;, ^ = CwA. (83) 

JM JM 

Integrating this by parts and using the relation a; = V x t;, we have 

Aa = - f (V X v ) d^cc + Intcj = - / p{b,v)(fx + Intg. (84) 

JM JM 

where Int^ denotes the integral over the surface S bounding M. This is the 
proposed Lagrangian A a in the physical space x, where 

pb = Vx* = Vx (CwA). (85) 

The Lagrangian of (83) is an integral of {CwA, uj) = CwA^[uj]. Given 
a vector field Y, the value taken on the vector Y measures the 

derivative of A^ (as one moves along the trajectory of W, i.e. the particle 
path) evaluated on a vector field Y frozen to the fiow generated by W. 

Using (82) in (83) and carrying out integration by parts, we obtain 

Aa= [ {A,EwH)d''x, (86) 

JM 

with -E'vk['^] defined in (76), where the integrated terms are dropped. Here, we 
redefine A a by this equation instead of (83). Hence, the expressions (83) and 
(84) must be modified by additional terms of surface integral over S, which 
are to be expressed symbolically with Intg. This Lagrangian (86) provides a 
source term of non-vanishing helicity. An example of such a fiow will be given 
in the last section 8.3. 

Note that, if f{t, x) is a scalar function, Cwf = dtf + Vf = Dtf- 
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7 Corrections to the previous Variation 



The Lagrangian of (84) should be included in the total Lagrangian to 
describe rotational flows of an ideal fluid. Accordingly, the Lagrangian density 
of the total Lagrangian (excluding surface terms) is deflned by 

Lt = LT[v,p,s,(j),ilj,b] = lp{v,v) - pe{p,s) - p {dt + v ■ 

-ps{dt + v-V)ip- p{b,v). (87) 

As before (Sec. 3.2), we take variations of the fleld variables f, p, s and 
potentials (f) and if). In addition, we must now include the variation of the new 
variable A. Independent variations are taken for those variables. Substituting 
the varied variables v + 5v,p + 5p, s + 5s, (j) + 6(f), ip + 6ip and A-\- SA into 
Lj'lv, p, s, (j), ip, A] and writing its variation as SL^^, we obtain 

6Lt = 5v ■ p{v -Vcj) - sVi) - w) - 5s p V)t^ (88) 
+ 5p {\u^ -h- Dt0 - s -V ■h) + {5 A, Ew^) 
+ 

where remaining terms denoted by are omitted since they are the same 

as (24). The new term w is deflned by 



-^5v' = -p{w,5v), (89) 

'w = {wi), w, = {d/dv'){v,h). (90) 

Thus, the variational principle, 5/ = for independent variations of 5v, 6p, 
6 A, etc. results in 

5v: V -V (p - sVi) - w = , (91) 

Sp: - h - Dtcp - sDtip - V ■ b = , (92) 

SA : EwOJ = , (93) 



where the omitted terms are the same as before. Thus, the continuity 

equation (29) and entropy equation (30) are obtained as before. The equation 
(91) gives an expression for the velocity v: 

V = V (f) + sVtp + w . (94) 

The vorticity a; = V x is represented by 

uj = Vs xV^p + V X w . (95) 
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The second term V x it? is a new term leading to non- vanishing vorticity in a 
homentropic fluid of uniform s. The equation (93) is the same as (77) by (76). 

Potential part and rotational part: 

The potential part of the present problem is described by (92) , and the rota- 
tional part by (77). In order to show that this is consisitent with the traditional 
formulation, suppose that the vorticity equation (93) is solved and a velocity 
field v{x, t) is found. In view of the relation u; = V x v, the velocity v should 
satisfy the following: 

dtV + u;xv = V H', (96) 
for a scalar function H'{x,t). Obviously, taking the curl of this equation re- 
duces to (77). 

On the other hand, according to Sec. 4. 2, we have the Euler's equation of mo- 
tion (52), because the new Lagrangian has no effect in the total Lagrangian 
of a space owing to its construction, and hence generates no new term. Using 
the identity (67), this reduces to the following: 

dtV + u;xv = -\/i^v^ + h). (97) 

Defining H' by H' = dt4> + s dtip + v ■ (b — w), the equation (92) is 

= - h - Dt(f) - sDti! - V ■ b = -h-dt(f) 

-s dtip-v- (V0 + sV^ + 6) = -h~H', (98) 

by using (94). In fact, it can be shown that both expressions (96) and (97) are 
equivalent with this H', as follows. First we note the following identities: 

{V xw) xv = {v- V)w + w'^Vv'' - V{v ■ w), (99) 
w'^Vv'' - {v ■ V) V0 = Vi^v^ - {v ■ V)0) - sdktpVv''. (100) 

The second is an identity with v given by (94). Using these, one can verify 
the following: 

dtW + (V X w) X V + Vlv ■ {w — b)] = dtV + {v ■ V)v + V/i. 

The right hand side vanishes by Eq.(52). Therefore, we obtain 

dtW + (y X w) XV = -V[v ■ {w - b)]. (101) 

It is not difficult to show that the equation (96) reduces to this equation for 
the above H' and the total velocity v of (94), where the total vorticity iv is 
given by (95). The equation (101) is a new equation to be satisfied by the 
rotational component w. 
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If s = So (uniform), then v = V$ + ($ = + So%jj), and oj = V x w. 
Subtracting (101) from (97), we obtain V[|^;^ + h + dt^ + v ■ {b - w)] = 0, 
which is equivalent to (98). 

Thus there is no contradiction between the traditional formulation and the 
present one. New aspects of the present formulation are described in the next 
section, in addition to the expression (101) for the rotational component w. 

8 Significance of A^i 

We consider here the significance of the Lagrangian Aa_ associated with rota- 
tion symmetry and its outcome derived from the variational principle. 

8. 1 Source of helicity 

Let us recall that pb is defined by (85). Using (82), the expression of ^ is 
rewritten as 

^ = dtA + v^dkA + AfcVi;*^ = dtA + {V X A)xv + V(A,.i;*^). 
Therefore, defining a vector B with S = V x A, we have 

pb = V x-^ = dtB + \/ x{Bxv) = Ew\B]. (102) 
Variation of of (84) with respect to a variation 5v is given by 

5Ka = - ( [(V X 5v) + (V X (5*, v)] d^x, 

JM 

where V x 5^ = V x [B x 5v) from (102) for the variation 5v on keeping A 
(hence B) fixed. Therefore, we have 

5Aa = - jj^j{(y X * + X b), 5v^d?x. 

This leads to the following expression for w defined by (89): 

pw = V + UJ X B = Ew[B]+uj X B, (103) 
where Ew[B] = dtB + V x {B xv) from (76). 

For flows of a homentropic fluid with constant entropy Sq, the velocity is given 
by (94): v = V ^ + w = (p + SQip). Therefore, the helicity is 

H = J uj-v d^x = j u-w d?x = j p"^uj- Ew[B] d^x, (104) 
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since V ■ cj = 0. Thus it is found that the term generates the hehcity. 

8.2 Uniqueness of transformation 

Transformation from the Lagrangian a space to Eulerian x{a) space is de- 
termined locally by nine components of the matrix dx^ /dal' . However, in the 
previous solution considered in Sec. 4. 2, we had three relations (48) between 
the cc- velocity v = {Xr,Yr, Zr) and the a- velocity {Va,Vb,Vc), and another 
three relations (51) between the a;-acceleration A = {X^-t, Y^-^, Z^-r) and the 
a-acceleration [Aa, Af,, A^). These six relations are not sufficient to determine 
the nine matrix elements dx^ / daK But now, we are equipped with a new struc- 
ture of rotational symmetry which is expected to resolve this indefiniteness. 

In fact, the remaining three conditions are found from the equation (72) 
connecting the a;-vorticity oj = {uJx,ujy,LJz) and the a-vorticity r2a(o.) = 
{fla,flb,^c)- From (72), these are 

Qa = {dbV dcZ - dcV dbz) + Uy {dbZ dcX - dcZ dbx) + uj^ (dbX dcV - dcX dby), 

(105) 

with other two equations for fib and flc determined cyclically. 

There are three vectors (velocity, acceleration and vorticity) determined by 
evolution equations subject to initial conditions in each space of x and a 
coordinates. Transformation relations of the three vectors suffice to determine 
the nine matrix elemets dx'^ /da} locally. Thus, the transformation between 
the Lagrangian a space and Eulerian x{a) space is determined uniquely. In 
this sense, the equation of the vorticity is essential for the uniqueness of the 
transformation between Lagrangian and Eulerian coordinates. 

8.3 Two linking vortices (an example) 

As an example of flows of non-zero helicity, we consider a particular flow 
field induced by tangling of two vortices, i.e. a rectilinear line- vortex L linked 
with a vortex ring R of infinitesimal vortex-cores, whose cross-sections are 
both circular and strengths are 71, and 7/? respectively. The helicity H is 
immediately found as (see below), hence it is non-zero. Therefore, this 

fiow field cannot be given a representation of a Clebsch-type form. 

We consider a fiow of an incompressible fiuid of density po in a cylindrical 
coordinate frame (x, r, 0). Suppose that we have a rectilinear line- vortex L 
of strength 7l(> 0) coinciding with the x-axis (axis of cylindrical symmetry) 
and a vortex ring R of strength 7i?(> 0) which coincides temporarily with a 
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circle of radius R in the plane x = centered at the origin. By symmetry 
consideration, it is readily seen for this initial configuration that the vortex 
ring R simply translates parallel to the x axis with a constant velocity U (say) 
without change of its form and its radius (apart from rotaional motion around 
the X axis without changing its form), while the rectilinear vortex L stays at 
the same position without change. With respect to this coordinate frame, the 
velocity field (at the initial instant) of the rectilinear vortex L is represented 
by a vector potential = {F{r), 0, 0), while that of the vortex ring R is given 
by \E'r = (0, 0, G{x, r)). Explicit expressions of velocity are Vl = V x l and 
Vji = 'V X "^ji respectively, and 



VL = iO,0,-drF), F = -^logr, 
Vn={r"'dr{rG),-d,G, 0), 



/o a/x^ + — 2rR cos (f) + R^ 
Their vorticities are 

u^L = (7L Ky)Kz). 0, 0), LJR = (0, 0, 7^ 6{r - R)6{x)), (106) 

where 6 is the delta function, and [y, z) denotes the cartesian coordinates in 
the plane perpendicular to the x-axis. Hence, the helicity is given by 



in 

By elementary calculus, we find that the integral of each term gives the same 
value "^lIr- Thus we obtain non-zero helicity H = 27^7/? (Moffatt 1969). 

The potential B = V x A can be determined in the following way. In reference 
to the coordinate frame moving with the vortex ring R, the velocity field is 
steady (and the vorticity too). In this frame, a uniform velocity (— [/, 0,0) is 
added, and total velocity is given by 

w = vl + vr+ {-U, 0,0) = V X x = v{x) = {u, V, w), (107) 

where x = {F{r) — ^Ur,0,G{x,r)) = xi^): independent of t. In the equation 
(103): PqV = E]y[B]+uj x B, the vorticity lj is nonzero only on the rectilinear 
line L and the circular curve R, and vanishes {uj = 0) at all other points. At 
points neither on L nor on R, we have the equation: 

Pov = Ew[B], Ew[B]=dtB + V X {B xv). (108) 

Let us seek a general solution to this equation for a given velocity field 
v{x). Suppose that we have a particle-path solution x = X{a,t) satisfy- 
ing dtX{a, t) = v{X), where Dta = for a = a{x, t) and a{x, 0) = x. Next, 
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we introduce a particular field B^{a,t) defined by 

which satisfies the equation: E]v[B^] = dtB^ + V x [B^ x v) = (Cauchy's 
solution) for an arbitrary initial field B^{x, 0). A general solution to the equa- 
tion (108) is given by 

B{x,t) = B,{a,t) + Boix), (109) 

where Bq{x) is a time-independent field and a particular solution to (108). In 
fact, substituting (109) into (108), we have poV x x = V x (Bq x v), in view 
of (107) and E\y[B^] = 0. Thus, we obtain 

Bo XV = \/i) + poX {= Q), 

where ip is a scalar function satisfying v ■ Q = {v ■ + poV • % = 0. Then, 
the solution Bq is given by 

Bo = (0, -Q^, Qr)/U, Q = {Qx, Qr, Q<j>)- 



9 Summary and discussion 

Following the scenario of the gauge principle of field theory, it is found that 
the variational principle of fluid motions can be reformulated successfully in 
terms of a covariant derivative and Lagrangians defined appropriately. The 
Lagrangians are determined such that a gauge invariance is satisfied under 
local Galilean transformations including the rotational symmetry. The fluid 
material is characterized thermodynamically in terms of mass density and 
entropy density. In addition, local rotation [i.e. vorticity) of the material is also 
taken into account by a newly-introduced Lagrangian K^- A gauge-covariant 
derivative Vf was defined in the present formulation including the rotational 
symmetry (it was up to §5 for the translation symmetry alone). Galilean 
invariance requires that the covariant derivative should be the convective time 
derivative following the motion of the background material, i.e. the so-called 
Lagrange derivative. The covariant derivative is an essential building block of 
the gauge theory. Applying Vt to the position coordinate x{a,t) of a fluid 
particle a, we obtain the particle velocity v. It is found that the vorticity 
uj is in fact a gauge field associated with the rotational symmetry. Lastly, 
an example of flow with non-vanishing helicity is presented together with a 
general solution of the vector potential for the rotational part. 

The variational principle {i.e. the action principle) in the space of particle 
coordinate a results in the Euler-Lagrange equation. This yields the Noether 
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thorem, i.e. conservation equations of energy-momentum. The equations thus 
obtained are the Lagrangian form of the equation of motion and energy equa- 
tion, which are transformed to Euler's equation of motion and the conservation 
equation of total energy ^v"^ + e in Eulerian space. 

The variational principle with respect to Eulerian space coordinate x results 
in (i) the expression of velocity v in terms of gauge potentials, (ii) an integral 
describing a potential part of the velocity, (iii) the equation of continuity, (iv) 
the entropy equation. In addition to these results obtained in the previous case 
of translation symmetry alone, new results in the present formulation are as 
follows, (v) The velocity v includes a new rotational term w, (vi) the vorticity 
equation (93) is deduced from the variational principle and (vii) the equation 
(101) is derived for the rotational component w. It is noteworthy that the new 
Lagrangian of (86) yields a source term for the helicity defined by (104). 

Transformation from the Lagrangian a space to Eulerian x is determined lo- 
cally by nine components of the matrix dx^/daK However, in the solution 
obtained previously (Kambe 2007) without taking account of rotational sym- 
metry, we had only six relations for the velocity (48) and acceleration (51) in- 
cluding the matrix element dx'^/daK Both of the equations are invariant under 
rotational transformations of a displacement vector AX = (AX, AY, AZ), so 
that the vector AX is not uniquely determined. In the improved formulation 
taking account of the rotational symmetry with three additional equations 
from (72), there are three vectors (velocity, acceleration and vorticity) deter- 
mined by evolution equations subject to initial conditions in each space of x 
and a. Transformation relations of the three vectors suffice to determine the 
nine matrix elements dx^/dd locally. Thus, the transformation between the 
Lagrangian a space and Eulerian x{a) space is determined uniquely. In this 
sense, the equation (105) for the vorticity is essential to the uniqueness of the 
transformation between Lagrangian and Eulerian coordinates. 

Accordingly, the present variational formulation is self-consistent and compre- 
hensively describes flows of an ideal fluid. 
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Appendix: Background of the theory 

This appendix reviews the background of the theory and describes the scenario 
of the gauge principle in a physical system. 



A Gauge invariances 

In the theory of electromagnetism, it is well-known that there is an invariance 
under a gauge transformation of electromagnetic potentials consisting of a 
scalar potential and a vector potential A. The electric field E and magnetic 
field B are represented as E = — V0 — {l/c){dA/dt), and B = V x A, where 
c is the light velocity. The idea is that the fields E and B are unchanged 
by the following transformation: (0, A) — (0', A'), where 0' = — c~^dtf, 
A' = A + V / with f{x,t) being an arbitrary differentiable scalar function of 
position X and time t. 

Gauge invariance is also known in certain rotational fiows of the Clebsch rep- 
resentation (Eckart 1960). Its velocity field is represented as (Sec. 3. 2): 

V = V(j) + sVilj, 

where 0, and s are scalar functions, and ip satisfies the equation Dftp = 0. 
Let 0', ip' and s' be a second set giving the same velocity field v, which implies 
the following: sVip — s'Vip' = V(0' — 0). This relation will hold true by the 
relation 0' — = ^(-0, if s = OF /dip and s' = —dF/dip'. Therefore, s and 
■0 are determined only up to such a contact transformation, and transforms 
by the addition of the generating function F{ip, It can be shown that the 
two sets of triplet give the identical fiow field. Thus, there exist some freedom 
in the expressions of physical fields in terms of potentials. 



B Related aspects in quantum mechanics and relativity theory 

In quantum mechanics, Schrodinger's equation for a charged particle of mass 
m and electromagnetic fields are invariant with respect to a gauge transforma- 
tion. This is as follows. In the absence of electromagnetic fields, Schrodinger's 
equation for a wave function tp of a free particle m is written as S[ip] = ih dtip — 
2^ Pfc = 0) where pk = —ihdk is the momentum operator {dk = d/dx^ for 
k = 1,2,3). In the presence of electromagnetic fields, Schrodinger's equation 
for a particle with an electric charge e is obtained by the transformation: 

dt^dt + {e/ih)Ao, dk^dk + {e/zhc)Ak, (A; = 1,2,3), (B.l) 
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where (A^) = {Aq, Ai, A2, A3) = {Aq, A) = (—0, A) is a four-vector poten- 
tial, actually a covector denoted with lower indices. A point in space-time 
is denoted by a four vector with upper indices, ( ) with 



X 



ct. Thus, we obtain the equation with electromagnetic fields: 



SaH'] = ihdti^ - e4)i: - — 

2m 



-th){dk + —Ak) ^ = 0. 

the J 



From the four potential 



A), we obtain E and B. 



Suppose that a wave function %Ij{x^) satisfies the equation Sa[ 
consider the following set of transformations of ip{x^) and A^: 



0, and 



xl)\x^) =exp[za(x^)] %Ij{x^), 

Afc ^ = + 80, ^ 0' = _ c-^dtf5, 

where a = {e/hc)(3. Then, it is shown readily that the transformed function 
il)'{x^) satisfies the Schrodinger equation S'a'['^'] = 0. This is the gauge invari- 
ance of the system of an electric charge in electromagnetic fields. The system 
is said to have a gauge symmetry. 

A relativistic Lagrangian [i.e. Lorentz- invariant) is defined for fluid flows by 
Soper (1976), and the energy-momentum tensor T^*^ is derived from it, where 
/i, 1/ = 0,1,2,3 with denoting the time component. In this relativistic for- 
mulation, the continuity equation and an energy equation of fluid flows are 
derived from the conservation of particle numbers and an equation of energy 
conservation d^T^" = 0. Expanding the energy density and energy flux 
T^'' with respect to an infinitesimal ratio (3 = v/c {k = 1,2, 3, where c is the 
light velocity and v a representative fluid velocity), the continuity equation is 
obtained from the lowest order terms (equivalent to the conservation of par- 
ticle number), while nonrelativistic energy conservation is obtained from the 
next order terms proportional to The entropy equation is another form of 
the energy conservation. 



C Brief scenario of gauge principle 



In the gauge theory of particle physics (Weinberg 1995; Frankel 1997), a free- 
particle Lagrangian Afree['0] is defined first for the wave function %l){x^) of a 
particle with an electric charge. Let us consider the following gauge transfor- 
mation: ip I— s> ip. If Afree is invariant under the transformation when a is 
a constant, it is said that Afree has a global gauge invariance. In spite of this, 
it often happens that Af^ee is not invariant for a function a = a{x), i.e. Afree 
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is not gauge-invariant locally. In this circumstance, it is instrumental to intro- 
duce a new field in order to acquire local gauge invariance. If the new field (a 
gauge field) is chosen appropriately, local gauge invariance can be recovered. 

In the section B, the /oca/ gauge-invariance was acquired by replacing 9^ with 

= 9^ + Af, 

(see (B.l)), where Af^ = {e/ihc)A^, and A^j,[x) is the electromagnetic potential 
(and termed a connection form in mathematics). The operator is called 
the covariant derivative. 

Thus, when the original Lagrangian is not locally gauge invariant, the principle 
of local gauge invariance requires a new gauge field to be introduced in order to 
acquire local gauge invariance, and the Lagrangian is to be altered by replacing 
the partial derivative with the covariant derivative including the gauge field. 
This is the Weyl's gauge principle. 

In mathematical terms, suppose that we have a group Q of transformations 
and an element g{x) E Q (for x E M with M a space where ip is defined) and 
that the wave function iplx) is transformed as ip' = g{x)ip. In the previous 
example, g{x) = e*"^"^^ and the group is ^ = U{1). Introducing the gauge field 
A allows us to define a covariant derivative D = d + A as a generalization of 
the partial derivative d that transforms as gD = g{d + A) = [d' + A')g. If we 
operates the right hand side on ip., we obtain {d' + A')g'il) = {d' + A') ip' = D'lp' 
where D' = d' + A'. Thus, we obtain D'lp' = g Dip, showing that Dip transforms 
in the same way as ip itself. 

In dynamical systems which evolve with the time t, such as the present case 
of fluid flows, replacement is to be made only for the t derivative: dt Dt = 
dt + Aix). 
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